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(54) Adaptive equalization and interpolated timing recovery in a sampled amplitude read channel 
for magnetic recording 

(57) A sampled amplitude read channel for mag- 
netic disk recording which asynchronously samples (24) 
the analog read signal, adaptively equalizes (BIOS) the 
resulting discrete time sample values according to a tar- 
get partial response, extracts synchronous sample val- 
ues through interpolated timing recovery (B100). and 
detects digital data from the synchronous sample val- 
ues using a Viterbi sequence detector (34) is disclosed. 
With interpolated timing recovery, the equalizer filter 
and its associated latency are removed from the timing 
recovery loop, thereby allowing a higher order discrete 
time filter and a lower order analog filter. 
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Description 

FIELD OF INVENTION 

5 The present invention relates to the control of magnetic disk storage systems for digital computers, particularly to 

a sampled amplitude read channel that employs asynchronous sampling of an analog read signal, adaptive discrete 
time equalization, and interpolated timing recovery 

CROSS REFERENCE TO RELATED APPLICATIONS AND PATENTS 

10 

This application is related to other co-pending U.S. patent applications, namely application serial numbers 
08/440,515 entitled "Sanrpted Amplitude Read Channel For Reading User Data and Embedded Servo Data From a 
Magnetic Medium," 08/341,251 entitled "Sampled Amplitude Read Channel Comprising Sample Estimation Equaliza- 
tion, Defect Scanning, Channel Quality, Digital Servo Demodulation. PID RIter for Timing Recovery, and DC Offset 

15 Control," 08/313,491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording," and 
08/533,797 entitled "Improved Fault Tolerant Sync Mark Detector For Sampled Amplitude Magnetic Recording." This 
application is also related to several U.S. patents, namely U.S. Pat. No. 5,359,631 entitled "Timing Recovery Circuit for 
Synchronous Waveform Sampling," 5,291,499 entitled "Method and Apparatus for Reduced-Complexity Viterbi-Type 
Sequence Detectors." 5,297.184 entitled "Gain Control Circuit for Synchronous waveform Sampling." 5.329,554 enti- 

20 tied "Digital Pulse Detector," and 5,424,881 entitled "Synchronous Read Channel." All of the above-named patent appli- 
cations and patents are assigned to the same entity, and all are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

25 In magnetic storage systems for computers, digital data serves to modulate the current in a read/write head coil in 
order to write a sequence of corresponding magnetic flux transitions onto the surface of a magnetic medium in concen- 
tric, radially spaced tracks at a predetermined baud rate. When reading this recorded data, the read/write head again 
passes over the magnetic medium and transduces the magnetic transitions into pulses in an analog read signal that 
alternate in polarity. These pulses are then decoded by read channel circuitry to reproduce the digital data. 

30 Decoding the pulses into a digital sequence can be performed by a simple peak detector in a conventional analog 
read channel or, as in more recent designs, by a discrete time sequence detector in a sampled amplitude read channel. 
Discrete time sequence detectors are preferred over simple analog pulse detectors because they compensate for inter- 
symbol interference (ISI) and are less susceptible to channel noise. As a result, discrete time sequence detectors 
increase the capacity and reliability of the storage system. 

35 There are several well known discrete time sequence detection methods including discrete time pulse detection 
(DPD). partial response (PR) with Viterbi detection, maximum likelihood sequence detection (MLSD). decision-feed- 
back equalization (DFE). enhanced decision-feedback equalization (EDFE). and fixed-delay tree-search with decision- 
feedback (FDTS/DF). 

In conventional peak detection schemes, analog circuitry, responsive to threshold crossing or derivative informa- 

40 tion, detects peaks in the continuous time analog signal generated by the read head. The analog read signal is "seg- 
mented" into bit cell periods and interpreted during these segments of time. The presence of a peak during the bit cell 
period is detected as a "1" bit, whereas the absence of a peak is detected as a "0" bit. The most common errors in 
detection occur when the bit cells are not correctly aligned with the analog pulse data. Timing recovery, then, adjusts 
the bit cell periods so that the peaks occur in the center of the bit cells on average in order to minimize detection errors. 

45 Since timing information is derived only when peaks are detected, the input data stream is normally run length limited 
(RLL) to limit the number of consecutive "0" bits. 

As the pulses are packed closer together on the concentric data tracks in the effort to increase data density, detec- 
tion errors can also occur due to intersymbol interference, a distortion in the read signal caused by closely spaced over- 
lapping pulses. This interference can cause a peak to shift out of its bit cell, or its magnitude to decrease, resulting in a 

so detection error. The ISI effect is reduced by decreasing the data density or by employing an encoding scheme that 
ensures a minimum number of "0" bits occur between "1" bits. For example, a (d.k) run length limited (RLL) code con- 
strains to d the minimum number of "0" bits between "1" bits, and to k the maximum number of consecutive "0" bits. A 
typical (1 .7) RLL 2/3 rate code encodes 8 bit data words into 12 bit codewords to satisfy the (1.7) constraint. 

Sanpled amplitude detection, such as partial response (PR) with Viterbi detection, allows for increased data den- 

55 sity by compensating for irrtersymbol interference and the effect of channel noise. Unlike conventional peak detection 
systems, sampled amplitude recording detects digital data by interpreting, at discrete time instances, the actual value 
of the pulse data. To this end. the read channel comprises a sampling device for sampling the analog read signal, and 
a timing recovery circuit for synchronizing the samples to the baud rate (code bit rate). Before sampling the pulses, a 
variable gain amplifier adjusts the read signal's amplitude to a nominal value, and a low pass analog filter filters the read 
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signal to attenuate aliasing noise. After sampling, a digital equalizer filter equalizes the sample values according to a 
desired partial response, and a discrete time sequence detector, such as a Viterbi detector, Interprets the equalized 
sample values in context to determine a most likely sequence for the digital data (i.e.. maximum likelihood sequence 
detection (MLSD)). MLSD takes into account the effect of ISI and channel noise in the detection algorithm, thereby 
5 decreasing the probability of a detection error. This increases the effective signal to noise ratio and. for a given (d.k) 
constraint, allows for significantly higher data density as compared to conventional analog peak detection read chan- 
nels. 

The application of sampled amplitude techniques to digital communication channels is well documented. See Y. 
Kabal and S. Pasupathy, "Partial Response Signaling". IEEE Trans. Common. Tech., Vol. COM-23. pp.921 -934. Sept. 

10 1 975; and Edward A. Lee and David G. Messerschmitt, "Digital Communication". Kluwer Academic Publishers. Boston, 
1990; and G.D. Forney, Jr., "The Viterbi Algorithm", Proc. IEEE, Vol. 61 . pp. 268-278, March 1973. 

Applying sampled amplitude techniques to magnetic storage systems is also well documented. See Roy D. Cideci- 
yan. Francois Dollvo. Walter Hirt, and Wolfgang Schott, "A PRML System for Digital Magnetic Recording", IEEE Journal 
on selected Areas in Communications. Vol. 10 No. 1, January 1992, pp.38-56; and Wood et al, "Vitertii Detection of 

15 Class IV Partial Response on a Magnetic Recording Channel IEEE Trans. Common., Vol. Com-34. No. 5. pp. 454- 
461. May 1986; and Coker Et al. "Implementation of PRML in a Rigid Disk Drive", IEEE Trans, on Magnetics, Vol. 27, 
No 6, Nov. 1991; and Carley et al, "Adaptive Continous-Time Equalization Followed By FDTS/DF Sequence Detec- 
tion". Digest of The Magnetic Recording Conference, August 15-17, 1994. pp. C3; and Moon et al. "Constrained-Com- 
plexity Equalizer Design for Rxed Delay Tree search with Decision Feedback". IEEE Trans, on Magnetics, Vol. 30, No. 

20 5. Sept. 1994; and Abbott et al, "Timing Recovery For Adaptive Decision Feedback Equalization of The Magnetic Stor- 
age Channel". Globecom'90 IEEE Global Telecommunications Conference 1990, San Diego, CA, Nov. 1990, pp.1794- 
1799; and Abbott et al, "Performance of Digital Magnetic Recording with Equalization and Offtrack Interference", IEEE 
Transactions on Magnetics, Vol. 27. No. 1 , Jan. 1991 ; and Cioffi et al, "Adaptive Equalization in Magnetic- Disk Storage 
Channels". IEEE Communication Magazine, Feb. 1990; and Roger Wood. "Enhanced Decision Feedback Equaliza- 

25 tion", fntermag'90. 

Similar to conventional peak detection systems, sampled amplitude detection requires timing recovery in order to 
correctly extract the digital sequence. Rather than process the continuous signal to align peaks to the center of bit cell 
periods as in peak detection systems, sampled amplitude systems synchronize the pulse samples to the baud rate. In 
conventional sampled amplitude read channels, timing recovery synchronizes a sampling clock by minimizing an error 
30 between the signal sample values and estimated sample values. A pulse detector or slicer determines the estimated 
sample values from the read signal samples. Even in the presence of ISI the sample values can be estimated and, 
together with the signal sample values, used to synchronize the sampling of the analog pulses in a decision-directed 
feedback system. 

A phase-locked-loop (PLL) normally implements the timing recovery decision-directed feedback system. The PLL 
35 comprises a phase detector for generating a phase error based on the difference between the estimated samples and 
the read signal samples. A PLL loop filter filters the phase error, and the filtered phase error operates to synchronize 
the channel samples to the baud rate. Conventionally, the phase en-or adjusts the frequency of a sampling clock which 
is typically the output of a variable frequency oscillator ( VFO). The output of the VFO controls a sampling device, such 
as an analog-to-digital (A/D) converter, to synchronize the sampling to the baud rate. 
40 As mentioned above, sampled amplitude read channels also commonly employ a discrete time equalizer filter to 
equalize the sample values into a desired partial response (PR4, EPR4. EEPR4. etc.) before sequence detection. To 
this end, adaptive algorithms have been applied to compensate in real time for parameter variations in the recording 
system and across the disk radius. For example. U.S. Pat. No. 5.381,359 entitied "Adaptation and Training of Digital 
Finite Impulse Response Filter Within PRML Sampling Data Detection Channel", discloses an adaptive equalizer filter 
45 that operates according to a well known least mean square (LMS) algorithm. 

Wk.i = W,-n*e,.X,. 

where W^ represents a vector of filter coefficients; is a programmable gain; ek represents a sample error between the 
50 filter's actual output and a desired output; and represents a vector of sample values from the filter input. In other 
words, the LMS adaptive equalizer filter is a closed loop feedback system that attempts to minimize the mean squared 
error between an actual output of the filter and a desired output by continuously adjusting the filter's coefficients to 
achieve an optimum frequency response. 

A problem associated with adaptive equalizer filters in sampled amplitude read channels is that the timing recovery 
55 and gain control loops can interfere with the adaptive feedback loop, thereby preventing the adaptive equalizer filter 
from converging to an optimal state. This non-convergence is manifested by the filter's phase and gain response drifting 
as it competes with the timing and gain control loops. An article by J.D. Coker et al. entitied "Implementation of PRML 
in a Rigid Disk Drive", published in IEEE Transactions on Magnetics, vol. 27. No. 6, Nov. 1991, suggests a three tap 
transversal filter comprising a fixed center tap and symmetric side taps in order to constrain the phase response of the 
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equalizer filter except in terms of a fixed group delay. Constraining the phase response of the adaptive equalizer in this 
manner, however, is a very sub-optimal method for attenuating interference from the timing recovery and gain control 
loops. Furthermore, it significantly reduces control over the adaptive filter's phase response, thereby placing the burden 
of phase compensation on the analog equalizer. 
5 Yet another problem associated with conventional adaptive equalizer filters is an inherent limitation on its order (i.e., ^ 

the number of coefficients): because the adaptive equalizer is inside the timing recovery feedback loop, its order must 
be limited to minimize the associated transport delay. Compensating for the deficiencies of the discrete time equalizer 
requires filtering the analog read signal with a higher order analog equalizer prior to the timing recovery loop, which is 
undesirable. 

10 There is, therefore, a need for an adaptive, discrete time equalizer filter in a sampled amplitude read channel hav- 
ing an improved method for constraining the phase and gain response in order to minimize interference from the timing 
recovery and gain control loops. A further aspect of the present invention is to remove the adaptive equalizer, and its 
associated latency, from the timing recovery loop, thereby allowing a higher order discrete time filter and a simplified 
analog filter. 

15 

SUMMARY OF THE INVENTION 

A sanpled amplitude read channel for magnetic disk recording which asynchronously samples the analog read sig- 
nal, adaptively equalizes the resulting discrete time sample values according to a target partial response, extracts syn- 

20 chronous sample values through interpolated timing recovery, and detects digital data from the synchronous sample 
values using a Viterbi sequence detector is disclosed. To minimize interference from the timing and gain control loops, 
the phase and magnitude response of the adaptive equalizer filter are constrained at a predetermined frequency using 
an optimal orthogonal projection operation as a modification to a least mean square (LMS) adaptation algorithm. Fur- 
ther, with interpolated timing recovery, the equalizer filter and its associated latency are removed from the timing recov- 

25 ery loop, thereby allowing a higher order discrete time filter and a lower order analog filter. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The above and other aspects and advantages of the present invention will be better understood by reading the fol- 
30 lowing detailed description of the invention in conjunction with the drawings, wherein: 

FIG. 1 is a block diagram of a conventional sampled amplitude recording channel. 

FIG. 2A shows an exemplary data format of a magnetic disk having a plurality of concentric tracks comprised of a 
plurality of user data sectors and embedded servo data sectors. 
35 FIG. 2B shows an exemplary format of a user data sector. 

FIG. 3 is a block diagram of the improved sampled amplitude read channel of the present invention comprising 
interpolated timing recovery for generating interpolated sample values and a synchronous data clock for clocking 
operation of a discrete time sequence detector and peak detector. 

FIG. 4A is a detailed block diagram of the prior art sampling timing recovery comprising a sampling VFO. 
40 FIG. 48 is a detailed block diagram of the interpolating timing recovery of the present invention comprising an inter- 
polator. 

FIG. 5 illustrates the channels samples in relation to the interpolated baud rate samples for the acquisition pream- 
ble- 

FIG. 6 shows an FIR fitter implementation for the timing recovery interpolator. 
45 FIG. 7 depicts a cost reduced implementation for the timing recovery interpolator. 

FIG. 8A is a block diagram of a conventional adaptive, discrete time equalizer filter in a sampled amplitude read 
channel- 

FIG. 8B shows the adaptive, discrete time equalizer of the present invention. 

FIG. 8C shows an alternative embodiment for the adaptive, discrete time equalizer of the present invention. 
50 FIG. 9A illustrates the present invention adaptive filter's gain response constrained at a normalized frequency of 
1/4T. 

FIG. 9B shovi^s the present invention adaptive filter's phase response constrained at a normalized frequency of 
1/4T 

FIG. 10 illustrates operation of an orthogonal projection operation of ttie present invention for constraining the gain 
55 and phase response of the adaptive filter. 

FIG. 1 1 A shows an implementation for a reduced cost orthogonal projection operation. 
FIG. 11 B shows an alternative embodiment for the reduced cost orthogonal projection operation. 
FIG. lie illustrates an innplementation for a gradient averaging circuit used in the reduced cost orthogonal projec- 
tion operation of FIG. 11 B. 
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FIG. 1 1 D shows yet another alternative embodiment for the reduced cost orthogonal projection operation of the 
present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

5 

Conventional Sampled Amplitude Read Channel 

Referring now to FIG. 1 , shown is a detailed block diagram of a conventional sampled amplitude read channel. Dur- 
ing a write operation, either user data 2 or preamble data from a data generator 4 (for example 2T preamble data) is 

10 written onto the media. An RLL encoder 6 encodes the user data 2 into a binary sequence b(n) 8 according to an RLL 
constraint. A precoder 1 0 precodes the binary sequence b(n) 8 in order to conrpensate for the transfer function of the 
recording channel 18 and equalizer filters to form a preceded sequence ~b(n) 12. The precoded sequence ~b(n) 12 is 
converted into symbols a(n) 16 by translating 14 ~b(N) = 0 into a(N) = -1. and --bCN) = 1 Into a(N) = +1 . Write circuitry 
9, responsive to the symbols a(n) 1 6, modulates the current in the recording head coil at the baud rate 1/T to record the 

15 binary sequence onto the media. A frequency synthesizer 52 provides a baud rate write clock 54 to the write circuitry 9 
and is adjusted by a channel data rate signal (CDR) 30 according to the zone the recording head is over. 

When reading the recorded binary sequence from the media, timing recovery 28 first locks to the write frequency 
by selecting, as the input to the read channel, the write clock 54 through a multiplexor 60. Once locked to the write fre- 
quency, the multiplexor 60 selects the signal 1 9 from the read head as the input to the read channel in order to acquire 

20 an acquisition preamble recorded on the disk prior to the recorded user data. A variable gain amplifier 22 adjusts the 
amplitude of the analog read signal 58, and an analog filter 20 provides initial equalization toward the desired response 
as well as attenuating aliasing noise. A sampling device 24 samples the analog read signal 62 from the analog filter 20, 
and a discrete time equalizer filter 26 provides further equalization of the sample values 25 toward the desired 
response. In partial response recording, for example, the desired response is often selected from Table 1. 

25 The discrete equalizer filter 26 may be implemented as a real-time adaptive filter which compensates for parameter 

variations over the disk radius (i.e.. zones), disk angle, and environmental conditions such as temperature drift. To this 
end. the filter 26 receives estimated sample values B143 generated by the timing recovery circuit 28: the estimated 
samples being input into an adaptive feedback loop and used to generate sample errors. The adaptive feedback loop 
conventionally employs a least mean square (LMS) algorithm to adapt the filter coefficients (i.e. . it adapts the frequency 

30 and phase response of the filter) until a minimum sample error is achieved. Operation of a conventional adaptive equal- 
izer filter is discussed in greater detail below. 

After equalization, the equalized sample values 32 are applied to a decision directed gain control 50 and timing 
recovery 28 circuit for adjusting the amplitude of the read signal 58 and the frequency and phase of the sampling device 
24, respectively. Timing recovery adjusts the frequency of sampling device 24 over line 23 in order to synchronize the 

35 equalized samples 32 to the baud rate. Frequency synthesizer 52 provides a course center frequency setting to the tim- 
ing recovery circuit 28 over line 64 in order to center the timing recovery frequency over temperature, voltage, and proc- 
ess variations. The channel data rate (CDR) 30 signal adjusts a frequency range of the synthesizer 52 according to the 
data rate for the current zone. Gain control 50 adjusts the gain of variable gain amplifier 22 over line 21 in order to match 
the magnitude of the channel's frequency response to the desired partial response. 

40 The equalized samples Y(n) 32 are also sent to a discrete time sequence detector 34, such as a maximum likeli- 
hood (ML) Viterbi sequence detector, which detects an estimated binary sequence ''b(n) 33 from the sample values. An 
RLL decoder 36 decodes the estimated binary sequence ''b(n) 33 from the sequence detector 34 into estimated user 
data 37. A data sync detector 66 detects the sync mark 70 (shown in FIG. 2B) in the data sector 15 in order to frame 
operation of the RLL decoder 36. In the atDsence of errors, the estimated binary sequence ''b(n) 33 matches the 

45 recorded binary sequence b(n) 8, and the decoded user data 37 matches the recorded user data 2. 

Data Format 

FIG. 2A shows an exemplary data format of a magnetic media comprising a series of concentric data tracks 13 
50 wherein each data track 1 3 comprises a plurality of sectors 1 5 with embedded servo wedges 1 7. A sen/o controller (not 
shown) processes the servo data in the servo wedges 17 and. in response thereto, positions the read/write head over 
a desired track. Additionally, the servo controller processes servo bursts within the servo wedges 17 to keep the head 
aligned over a centerline of the desired track while writing and reading data. The servo wedges 17 may be detected by 
a simple discrete time pulse detector or by the discrete time sequence detector 34. If the sequence detector 34 detects 
55 the servo data, then the format of the servo wedges 17 includes a preamble and a sync mark, similar to the user data 
sectors 15. 

FIG. 2B shows the format of a user data sector 15 comprising an acquisition preamble 68, a sync mark 70. and 
user data 72. Timing recovery uses the acquisition preamble 68 to acquire the correct sampling frequency and phase 
before reading the user data 72, and the sync mark 70 demarks the beginning of the user data 72 (see co-pending U.S. 
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patent application 08/313.491 entitled "Improved Timing Recovery For Synchronous Partial Response Recording"). 

To increase the overall storage density, the disk is partitioned into an outer zone 11 comprising fourteen data sec- 
tors per track, and an inner zone 27 comprising seven data sectors per track. In practice, the disk is actually partitioned 
into several zones with a different number of sectors in each zone, and the data recorded and detected at a different 
5 data rate in each zone. 

Improved Sampled Amplitude Read Channel 

FIG. 3 shows the improved sampled arrplitude read channel of the present invention wherein the conventional 
10 sampled timing recovery 28 of FIG. 1 has been replaced by interpolated timing recovery B100. In addition, the write 
frequency synthesizer 52 generates a baud rate write clock 54 applied to the write circuitry 9. and an asynchronous 
read clock 54 for docking the sampling device 24, the discrete time equalizer filter 26, and the interpolated timing recov- 
ery B100 at a frequency relative to the current zone (CDR 30). In an alternative embodiment, a first frequency synthe- 
sizer generates the write clock, and a second frequency synthesizer generates the read clock. 
15 The discrete equalizer filter B103 is real-time adaptive, receiving interpolated sample values 8102 and estimated 
sample values B143 from the interpolated timing recovery circuit B100 for use in a modified least mean square (LMS) 
algorithm which constrains the filter's gain and phase response according to the present invention, the details of which 
are set forth below. 

The interpolated timing recovery B100 interpolates the equalized sample values 32 to generate interpolated sam- 
po pie values B102 substantially synchronized to the data rate of the current zone. A discrete time sequence detector 34 
detects an estimated binary sequence 33 representing the user data from the interpolated sample values B102. The 
interpolated timing recovery B100 circuit generates a synchronous data clock B104 for clocking operation of the gain 
control 50. discrete time sequence detector 34, sync mark detector 66 and RLL decoder 36. 

25 Conventional Timing Recovery 

An overview of the conventional sampling timing recovery 28 of FIG. 1 is shown in FIG. 4 A, The output 23 of a var- 
iable frequency oscillator ( VFO) B1 64 controls the sampling clock of a sampling device 24 which is typically an analog- 
tCKiigital converter (A/D) in digital read channels. A multiplexor B159 selects the unequalized sample values 25 during 

30 acquisition and the equalized sample values 32 during tracking, thereby removing the discrete equalizer filter 26 from 
the timing loop during acquisition in order to avoid its associated latency. A phase error detector B1 55 generates a 
phase error in response to the sample values received over line B149 and estimated sample values -Yj^ from a sample 
value estimator B141, such as a slicer in a d=:0 PR4 read channel, over line B143. A loop filter B160 filters the phase 
error to generate a frequency offset B1 67 that settles to a value proportional to a frequency difference between the 

35 sampling clock 23 and the baud rate. The frequency offset A/ B167, together with the center frequency control signal 
64 from the frequency synthesizer 52, adjust the sampling clock 23 at the output of the VFO B164 in order to synchro- 
nize the sampling to the baud rate. 

A zero phase start B162 circuit suspends operation of the VFO B164 at the beginning of acquisition in order to min- 
imize the initial phase error between the sampling clock 23 and the read signal 62. This is achieved by disatjiing the 

40 VFO B164, detecting a zero crossing in the analog read signal 62, and re-enabling the VF0 164 after a predetermined 
delay between the detected zero crossing and the first baud rate sample. 

The estimated sample values B143 at The output of the slicer B141 are also input into the discrete time equalizer 
fitter 26 of FIG. 1 for use in a conventional least mean square (LMS) adaptation algorithm as is described in more detail 
below. 

45 

Interpolated Timing Recovery 

The interpolated timing recovery B100 of the present invention is shown in FIG. 48. The VFO B164 in the conven- 
tional timing recovery of FIG. 4 A is replaced with a modulo-Ts accumulator B120 and an interpolator B122. In addition. 

50 an expected sample value generator B151, responsive to interpolated sample values B102. generates expected sam- 
ples %^ used by the phase error detector B155 to compute the phase error during acquisition. A multiplexor B153 
selects the estimated sample values -Y^^.^ from the slicer B141 for use by the phase error detector 81 55 during track- 
ing. The data clock B104 is generated at the output of an AND gate B126 in response to the sampling clock 54 and a 
mask signal B1 24 from the modulo-Ts accumulator B1 20 as discussed in further detail below. The phase error detector 

55 B155 and the slicer 8141 process interpolated sample values B102 at the output of the interpolator B122 rather than 
the channel sample values 32 at the output of the discrete equalizer fitter 26 as in FIG. 4A. A PID loop filter B161 con- 
trols the closed loop frequency response similar to the loop filter B160 of FIG. 4A. 

The interpolated sample values Yk+^ B1Q2 and the estimated sample values -Yk+t from the slicer B141 are input 
into the adaptive, discrete equalizer fitter B103 of FIG. 3 for use by a modified least mean square (LMS) algorithm, the 
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details of which are set forth below. 

In the interpolated timing recovery of the present invention, locking a VFO to a reference frequency before acquir- 
ing the preannble is no longer necessary; multiplexing 60 the write dock 54 into the analog receive filter 20 (as in FIG. 
1) is not necessary Further, the sampling device 24 and the disaete equalizer filter 26, together with their associated 

5 delays, have been removed from the timing recovery loop; it is not necessary to multiplex B159 around the equalizer 
filter 26 between acquisition and tracking. However, it is still necessary to acquire a preamble 68 before tracking the 
user data 72. To this end. a zero phase start circuit B163 minimizes the initial phase error between the interpolated 
sample values and the baud rate at the beginning of acquisition similar to the zero phase start circuit B162 of FIG. 4A. 
However, rather than suspend operation of a sampling VFO 81 64. the zero phase start circuit 81 63 for interpolated tim- 

10 ing recovery computes an initial phase error t from the A/D 24 sample values 25 and loads this initial phase error into 
the modulo-Ts accumulator 8120. 

For more details concerning the RID loop filter 81 61 . phase en-or detector 81 55. expected sample generator 81 51 , 
and slicer 8141 . refer to the above referenced co-pending U.S. patent applications "Sampled Amplitude Read Channe|l 
Comprising Sample Estimation Equalization. Defect Scanning, Channel Quality. Digital Servo Demodulation, RID Filter 

15 for Timing Recovery, and DC Offset Control" and "Improved Timing Recovery For Synchronous Partial Response 
Recording • A detailed description of the modulo-Ts accumulator 8120, data clock 8104. and interpolator 8122 is pro- 
vided in the following discussion. 

Interpolator 

20 

The interpolator 8122 of FIG. 4B is understood with reference to FIG. 5 which shows a sarrpied 2T acquisition pre- 
amble Signal 8200 The target synchronous sample values 8102 are shown as black circles and the asynchronous 
channel sample values 32 as vertical arrows. Beneath the sampled preamble signal is a timing diagram depicting the 
corresponding timing signals for the sampling clock 54. the data clock 81 04 and the mask signal 81 24. As can be seen 
25 in FIG 5. the pi eamble signal 8200 is sampled slightly faster than the baud rate (the rate of the target values). 

The function of the interpolator is to estimate the target sample value by interpolating the channel sample values. 
For iliustratve purposes, consider a sinnple estimation algorithm, linear interpolation: 

Y(N-1) = x(N-1) + X • (x(N) - x(N-1)); where: (1) 

30 

x(N-i) and x(N) are the channel samples surrounding the target sample; and x is an interpolation interval proportional 
to a time drtfei ence between the channel sample value x(N-1) and tiie target sample value. The interpolation interval t 
is generated at the output of modulo-Ts accumulator 8120 which accumulates the frequency offset signal Af 8167 at 
the output ol the PID kx)p filter 8161 : 

35 

X = (ZAO MOD TS; where: (2) 

Ts IS the samp*ing period of the sampling clock 54. Since tiie sanrpling clock 54 samples the analog read signal 62 
slightly tastei irvin the tiaud rate, it is necessary to mask the data clock every time the accumulated frequency offset Af, 

40 integer divided b> Ts. increments by 1. Operation of the data clock 8104 and the mask signal 8124 generated by the 
modulo-Ts accunvitator 8120 is understood with reference to the timing diagram of FIG. 5. 

Assigning me mterpolator implements the simple linear equation (1) atxrve, then channel sample values 8202 and 
8204 aie us^sO to gtfrnsiate the interpolated sample value corresponding to target sample value 8206. The interpolation 
interval i 8208 is geoeialed according to equation (2) above. The next interpolated sample value corresponding to the 

45 next target yatue 821 0 is computed from channel sample values 8204 and 821 2. This process continues until the inter- 
polation intef va! i 8214 would be greater than Ts except that it "wraps" around and is actually t 8216 (i.e., the accumu- 
lated frequency onset A', integer divided by Ts. increments by 1 causing the mask signal 81 24 to activate). At this point, 
the data ck)cK Bl04 is masked by mask signal B124 so that the interpolated sample value corresponding to the target 
sample value B220 is computed from channel sample values 8222 and 8224 rather than channel sample values 8218 

50 and 8222 

The simple linear interpolation of equation (1) will only work if the analog read signal is sampled at a much higher 
frequency than the baud rate. This is not desirable since operating the channel at higher frequencies increases its com- 
plexity and cost Therefore, in the preferred embodiment the interpolator 8122 is implemented as a filter responsive to 
more than two chiannel samples to compute the interpolated sample value. 
55 The ideal disaete time phase interpolation filter has a flat magnitude response and a constant group delay of x: 

C,(eh = ej"^ (3) 

which has an ideal impulse response: 
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sinc(7i • (n-x/Tg)). (4) 

Unfortunately, the above non-causal infinite impulse response (4) cannot be realized. Therefore, the impulse response 
of the interpolation filter is designed to be a best fit approximation of the ideal impulse response (4). This can be accom- 
5 plished by minimizing a mean squared error between the frequency response of the actual interpolation filter and the 
frequency response of the ideal interpolation filter (3). This approximation can be improved by taking into account the 
spectrum of the input signal, that is. by minimizing the mean squared error between the input spectrum multiplied by 
the actual interpolation spectrum and the input spectrum multiplied by the ideal interpolation spectrum: 

10 C,(e*'")X(e^") - C,(e^")X(e^"); where: (5) 

Ci.(e^"*) is the spectrum of the actual interpolation filter; and X(e'^) is the spectrum of the input signal. From equation 
(5), the mean squared error is represented by:^ 



15 



:^ = |C,(e^^'^)- e^^'^^|'|X(e^")|'da,; where: (6) 



20 X{e^^) is the spectrum of the read channel (e.g., PR4, EPR4, EEPR4 of Table 1 or some other partial response spec- 
trum). 

In practice, the above mean squared error equation (6) is modified by specifying that the spectrum of the input sig- 
nal is bandlimitted to some predetermined constant 0 < co < aTc where 0 < a < 1 ; that is: 

25 |X(e*")| = 0, forH^ ct^ 

Then equation (6) can be expressed as: 



30 



= £ J IC.(e'")- e''-"|'|X(eni^da,. (7) 



The solution to the minimization problem of equation (7) involves expressing the actual interpolation filter in terms of its 
35 coefficients and then solving for the coefficients that minimize the error in a classical mean-square sense. 
The actual interpolation filter can be expressed as the FIR polynomial: 

n=R-1 

5^(e^'')= 52 ^t(n)e"^"": where: (8) 

40 n=-R 

2R is the number of taps in each interpolation filter and the sample period Ts has been normalized to 1 . A mathematical 
derivation for an interpolation filter having an even number of coefficients is provided below. It is within the ability of 
45 those skilled in the art to modify the mathematics to derive an interpolation filter having an odd number of coefficients. 
Substituting equation (8) into equation (7) leads to the desired expression in terms of the coefficients C^(n): 

= i T l"i;'c,(n)e-^-" - e^-|^|X(ehl'c^a>. (9) 

The next step is to take the derivatives of equation (9) with respect to the coefficients C^(n) and set them to zero: 

= Oforno = -R 0,1 R-1. (10) 



^C,(n^) 

After careful manipulation, equation (10) leads to: 
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J E <^T(n)cos(a>(n^-n)) .cos(<o(n^+x)) |X(e^"*)|2c/(o = 0 



(11) 



forno = -R, 



Defining ^{r) as: 



0. 1 R-1 



10 



<t»(r)= J |X(e'"*)r cos(a)r) d(o 



(12) 



25 



3D 



and substituting equation (12) into equation (11) gives: 



n=R-i 

E C,(n)4»(n-n^) = ,|>(n^ + T) 



(13) 



for no = -R, .... 0, 1, R-1 . 

Equation (13) defines a set of 2R linear equations in terms of the coefficients C^{n). Equation (13) can be expressed 
more compactly in matrix form: 



is a column vector of the form: 



c^T is a Toeplitz matrix of the form: 



<t>jC^ = <D^; where: 



C, = [c,{-R), ....c,(0) c,(R.1)]' 



35 



40 



O-r = 



<t)(0) 
<t){l) 

<|>(2R-1) 



<1>(1) 
(l>(0) 



<}»(2R-1 



<I>(0) 



45 



and <I>^ is a column vector of the form: 



<l>^ = [<t>(-R + x). .... <t>(t). 4»(1 + x). .... (t>(R-1 + t)] *. 



The solution to equation (14) is: 



(14) 



= <I> J * a>^; where: 



(15) 



<l>j ^ is an inverse matrix that can be solved using well known methods. 



Table B2 shows example coefficients C^{n) calculated from equation (15) with 2R=6. a=0.8 and X(e^"')=PR4 . The 
implementation of the six tap FIR filter is shown in FIG. 6. A shift register B250 receives the channel samples 32 at the 
sampling clock rate 54. The filter coefficients C^(n) are stored in a coefficient register file B252 and applied to corre- 
55 sponding multipliers according to the current value of x B128. The coefficients are multiplied by the channel samples 
32 stored in the shift register B250. The resulting products are summed B254 and the sum stored in a delay register 
8256. The coefficient register file 8252 and the delay register 8256 are clocked by the data clock 8104 to implement 
the masking function described above. 

In an alternative embodiment not shown, a plurality of static FIR filters, having coefficients that correspond to the 
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different values of x, filter the sample values in the shift register B250. Each filter outputs an interpolation value, and the 
current value of the interpolation interval x B128 selects the output of the corresponding filter as the output B102 of the 
interpolator B122. Since the coefficients of one filter are not constantly updated as in FIG. 6, this multiple filter embod- 
iment increases the speed of the interpolator 8122 and the overall throughput of the read channel. 

5 

Cost Reduced Interpolator 

Rather than store all of the coefficients of the interpolation filters In memory, in a more efficient, cost reduced imple- 
mentation the coefficient register file 8252 of FIG. 6 computes the filter coefficients C^(n) in real time as a function of t. 
10 For example, the filter coefficients C^(n) can be computed In real time according to a predetermined polynomial in x 
(see, for example, U.S. Pat. No. 4,866.647 issued to Farrow entitled. "A Continuously Variable Digital Delay Circuit," the 
disclosure of which is hereby incorporated by reference). An alternative, preferred embodiment for computing the filter 
coefficients in real time estimates the filter coefficients according to a reduced rank matrix representation of the coeffi- 
cients. 

75 The bank of f iltei coefficients stored in the coefficient register file 8252 can be represented as an MxN matrix Aj^xN» 
where N is the depth of the interpolation filter (i.e., the number of coefficients C^(n) in the impulse response computed 
according to equation (15)) and M is the number of interpolation intervals (i.e., the number of x Intervals). Rather than 
store the entire matrix in memory, a more efficient, cost reduced implementation is attained through factorization 
and singular value decomposition (SVD) of the A^xn matrix. 

20 Consider that the Aj^^n matrix can be factored into an F|^xN ^NxN rnatrix. 

^ MxN = ^ MXN * ^ NxN * 

Then a reduced rank approximation of the A^^^N matrix can be formed by reducing the size of the Fj^j^n and Gnxn matri- 
25 ces by replacing N with L where L<N and. preferably. L«N. Stated differently, find the F^^l and C^jdsi matrices whose 
product best approximates the A^^^n nf^atrix. 

A MxN ~ ^ MxL * ^ LxN - 

30 The convolution process of the interpolation filter can then be carried out. as shown in FIG 7. by implementing the Glxn 
matrix as a bank of FIR filters 8260 connected to receive the channel sample values 32, and the F^^l matrix imple- 
mented as a kx>ki43 table 8262 indexed by x 8128 (as wit! become more apparent in the following discussion). Those 
skilled in the ait will lecognize that, in an alternative embodiment, the A^xu if^atrix can be factored into more than two 
matrices (i e . A = FGH ) 

35 The preferred metfxxJ tor finding the f^^i and Gy^^ matrices is to minimize the following sum of squared enrors: 

M n 

MxL * ^LxN)jn) 

(16) 

j= 1n= 1 

40 

The solution to equaton ( 16) can be derived through a singular value decomposition of the A^xn matrix, comprising the 
steps of 

45 1 . performing an GVD on the Af^^^ matrix which gives the following unique factorization (assuming M>N): 

A MxN = U Mxt4 • D NXN • V^xN where: 

Um.n * ^ unitary matrix; _ . 

50 D^.N a N)iN diagonal nr^trix {a^ . .. ctn) where a; are the singular values of A|^xN' arKi ^ aa - ^ cfn ^ 

0: and 

^NxN ^ a NxN unitary nr^trix; 

2. selecting a predetermined L number of the largest singular values a to generate a reduced size diagonal matrix 

55 Dlxl: 
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0 .. 
0 0.. 



0 



0 ... 0 Or 



3. extracting the first L columns from the U^xn matrix to form a reduced U^xl matrix 



25 



u 



MxL - 



u 



1.1 



u 



M, 1 



u 



1,L 



u 



M, L 



U 



l.N 



... u 



M.N 



4. extracting the first L rows from the Vn^n matrix to form a reduced Vlxn matrix: 



55 



40 



V, 



LxN - 




5. defining the F^xl and Glxn matrices such that: 



BO (for example, let F ^xl = U mxl * D and G = V lxn )• 

In the above cost reduced polynomial and reduced rank matrix embodiments, the interpolation filter coefficients 
C,(n) are computed in real time as a function of t; that is. the filter's impulse response h(n) is approximated according to: 



h(n.t) = c ,{n) = ^ G ((n) • f(i.t) : where: 



(17) 
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f(i.x) is a predetermined function in x (e.g., polynomial in t or x indexes the above ff^xi Tiafrix); L is a degree which 
determines the accuracy of the approximation (e.g., the order of the polynomial or the column size of the above F^^l 
matrix); and Gi(n) is a predetermined matrix (e.g.. the coefficients of the polynomial or the above Glxn matrix). As L 
increases, the approximated filter coefficients C^(n) of equation (17) tend toward the ideal coefficients derived from 
5 equation (1 5). It follows from equation (1 7) that the output of the interpolation filter Y(x) can be represented as: 

N L 

Y(x)= 2^ U(x-n) ^Gj(n) • f(i.x) (18) 

n= 1 i= 1 

10 

where U(x) are the channel sample values 32 and N is the number of interpolation filter coefficients C^(n). 

Referring again to FIG. 6. the coefficient register file can compute the interpolation fitter coefficients C^{n) according 
to equation (17) and then convolve the coefficients CT;(n) with the channel samples U(x) 32 to generate the interpolated 
sample values B102 synchronized to the baud rate. However, a more efficient implementation of the interpolation filter 
15 can be achieved by rearranging equation (18): 

L N 

Y(x)= "Z^Clx) • U(x-n) (19) 

i = 1 n = 1 

20 

FIG. 7 shows the preferred embodiment of the interpolation filter according to equation (19). In the polynomial 
embodiment, the function of x is a polynomial in x. and the matrix Gi(n) are the coefficients of the polynomial. And in the 
reduced rank matrix embodiment, the function of x Is to index the above Fi^^l nnatrix B262, and the second summation 
25 in equation (19). 

N 

^ Gi(n) . U(x-n) 
n= 1 

30 

is implemented as a bank of FIR filters B260 as shown in FIG. 7. Again, in equation (19) L is the depth of the approxi- 
mation function f(i.x) (e.g.. the order of the polynomial or the column size of the above Ffji^L matrix) and N is the depth 
of the interpolation filter's impulse response (i.e., the number of coefficients in the impulse response). It has been deter- 
35 mined that N=8 and L=3 provides the best performance/cost balance; however, these values may increase as IC tech- 
nology progresses and the cost per gate decreases. 

Conventional Adaptive Equalizer 

40 FIG. 8A illustrates a prior art adaptive, discrete time equalizer that operates according to the well known least mean 
square (LMS) algorithm. 

Wk.i =W,-M-eK-X,. 

45 or alternatively, 

where represents a vector of FIR filter coefficients; >i is a programmable gain; e^ represents a sample error (or vec- 
50 tor of sample errors e^J between the FIR filter's actual output and a desired output; and X,^ represents a vector of sam- 
ples values (or a scalar XjJ from the FIR filter input. To better understand operation of the present invention, the secor^ 
representation of the LMS algorithm is used throughout this disclosure. 

The desired filter output is the estimated sample values - at the output of slicer B141 . The estinnated sample val- 
ues -Yk are subtracted from the FIR filter's output Y^ to generate the sample error ek- The LMS algorithm attenrpts to 
55 minimize the sannple error in a least mean square sense by adapting the FIR filter coefficients; that is. it adjusts the FIR 
filter's gain and phase response so that the overall channel response adequately tracks the desired partial response 
(e.g., PR4, EPR4, EEPR4. etc.). 

As previously mentioned, interference from the timing recovery 28 and gain control 50 loops can prevent the adap- 
tive, discrete time equalizer 26 from conve-ging to an optimal state. For example, a phase adjustment in the adaptive 
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filter 26 can affect the sampling phase error for timing recovery 28. Timing recovery 28 compensates for the filter's 
phase adjustment by adjusting its sampling phase; this adjustment can result in yet another phase adjustment by the 
adaptive equalizer 26. Thus, the phase response of the adaptive equalizer may never converge. Similarly the gain con- 
trol loop 50 can interfere with the gain response of the adaptive filter 26 and prevent it from converging. 

5 

Constrained Adaptive Equalizer 

FIG. SB illustrates operation of the adaptive, discrete time equalizer filter B1 03 of the present invention. Sample val- 
ues from the A/D 24 are input over line 25 to a discrete time FIR filter CI 00 comprising a predetermined number of coef- 
10 _ ficients. the values of which determine the filter's gain and phase response. Because the FIR filter C100 operates on 
the sample values prior to the interpolated timing recovery loop B1 00, its order can be increased over the prior art with- 
out adversely affecting the latency of timing recovery (i.e., the number of filter coefficients can be increased). 

The output Yk 32 of the FIR filter CI 00 is input into the interpolator 8122 for generating the interpolated sample 
values Yk^.^ B102. The interpolated sample values Y^^^ 8102 are input into a slicer 8141 (FIG. 4B) which generates 
75 estimated sample values ~Yk+^. The estimated sample values -Vy^^^ are subtracted from the interpolated sample val- 
ues Yk+^ at adder C102 to generate a sample error value e^^^ C104 that is synchronized to the baud rate rather than 
the sample rate. Because the LMS algorithm operates on sample values at the sanrple rate, it is necessary to con- 
vert the error value ek+t CI 04 into an error value e^ Cl 1 2 synchronous to the sample rate. This is accomplished by an 
interpolation circuit 0106 which computes an interpolated error value e^ C112 from the baud rate error values Bi^^^ 
20 CI 04. Preferably the eror value interpolation circuit Cl 06 is implemented as a first order linear interpolation, but it may 
be a simple zero order hold, or a more complex interpolation filter as described above. 

FIG. 80 shows an alternative embodiment for generating the error value eK- As illustrated, the estimated ideal sam- 
ple values ~ Yk+^ from the slicer 81 41 are interpolated by the interpolation circuit 01 06 to generate estimated ideal sam- 
ple values --Yk CI 20 which are subtracted 0102 from the equalized sample values Yk 32 at the output of the FIR filter 
25 01 00 to generate the error value e^ 

In both embodiments, the error value e^ 0112 is input into a modified LMS circuit 0114 which computes updated 
filter coefficients Wk 01 1 6 according to. 

Wk.i - Wk-M-Pv,V2-'-(X,.ek) 

30 

where the operation PviV2-'- is an orthogonal projection operation which constrains the gain and phase response of the 
FIR filter 01 00 in order to attenuate interference from the gain and timing loops. 

Operation of the orthogonal projection operation Pv-,V2^ will now be described in relation to the gain and phase 
response of the FIR filter 0100. FIG.s 9A and 9B show an example gain and phase response, respectively for the FIR 
35 filter 01 00 of the present invention. The gain and phase response vary over time as the filter adapts to parameter 
changes in the recording system: that is, the filter continuously adjusts the channePs overall frequency response so that 
it matches the desired partial response (PR4. EPR4. EEPR4, etc.) as best possitrfe. In the present invention, interfer- 
ence from the timing recovery and gain control loops is attenuated by constraining the gain and phase response of the 
FIR filter C100 at a predetermined frequency using an orthogonal projection operation FV'^V2^. 
40 Referring again to FIG 9A, the gain (magnitude) response of the FIR filter 0100 has been constrained to a prede- 
termined value (denoted by g) at the normalized frequency of 0.5 (l/4Ts). Similarly, the phase response of the FIR filter 
OlOO has been constrained to kn: at the normalized frequency of 0.5 as shown in FIG. 9B. In effect, these constraints 
allow the gain and phase to vary (adapt) at all frequencies except at the normalized frequency of 0.5, thereby constrain- 
ing the filter's frequency response in a manner that attenuates interference from the gain and timing loops. The gain 
45 constraint g is relatively arbitrary except that it is selected to optimize the dynamic range of the filter's coefficients. How- 
ever, constraining tiie filter's response at the normalized frequency of 0.5 and selecting a phase constraint of Ictc 
reduces the complexity of the ortiiogonal projection operation Pv^ V2-^, and simplifies implementation of the zero phase 
start circuit 8163 (FIG. 4B) and sync mark detector 66 (FIG. 1). 

As mentioned above, the zero phase start circuit 81 63 (FIG. 4B) minimizes the initial phase error between the inter- 
so polated sample values and the baud rate at the beginning of acquisition by computing an initial phase error x from the 
A/D 24 sample values 25 and then loading this initial phase error into the moduIo-Ts accumulator 8120. To compute 
the initial phase error x. the z&ro phase start circuit 81 63 must take into account the phase delay of the adaptive equal- 
izer filter 81 03 since interpolated timing recovery 81 00 operates on the equalized samples 32, not the A/D 24 samples 
25. With the acquisition preamble 68 (FIG. 2B) having a frequency of 1/4T (i.e.. 0.5 normalized), constraining the phase 
55 response of the adaptive equalizer 8103 at the preamble frequency (1/4Ts) fixes the phase delay of the equalizer 8103 
during acquisition, thereby allowing the zero phase start circuit B163 to accurately compute the initial phase error t. 
Furthermore, since the phase constraint at the preamble frequency is fixed at Ktc, the phase delay of the equalizer 8103 
will either be zero or 180° (i.e., the adjustment to the initial phase error is nothing or a sign change). 

Constraining the phase response of the adaptive equalizer 81 03 to kn at the preamble frequency also simplifies 



13 



WSDOCID: <EP_._080&447A2_L> 



EP 0 805 447 A2 



10 



implementation of the sync mark detector 66 (FIG. 1) in sampled amplitude read channels that use the conventional 
synchronous sanrpling timing recovery 28 (FIG. 4A) rather than interpolated timing recovery B100 (FIG. 4B). Operation 
of the sync mark detector 66 is described in detail in the above referenced co-pending U.S. Patent Application 
08/533.797 entitled "Improved Fault Tolerant Sync Mark Detector For sanrpled Amplitude Magnetic Recording." As 
described therein, the sync mark detector 66 is enabled coherent with the end of the acquisition preamble 68 and rela- 
tive to the transport delay from the output of the A/D 24 to the sync mark detector 66. With an adaptive equalizer, the 
transport delay will vary unless the filter's phase response is constrained at the acquisition preamble frequency by. for 
exanple, using an orthogonal projection operation PviV2-^ of the present invention. 

Turning now to the implementation details of the orthogonal projection operation Pv-| V2-^, the equalizer's frequency 
response is 



C(e''"*) = ZCKe 



-jk2nfT 



where are the coefficients of the equalizer's impulse response. At the preamble frequency (1/4T). the equalizer's fre- 
quency response is 



25 where the sampling period has been normalized to T=1 . In matrix form, the equalizer's frequency response at the pre- 
amble frequency is. 



30 



35 



40 



(•")■■ 



(•'■) 



(j)' 
(j) 



-1 



(j) 





1 




-j 
-1 




j 


1 





Those skilled in the art will recognize that shifting the time base will lead to four different, but functionally equivalent, 
frequency responses at the preamble frequency (i.e., [1.-j.-1 j,- -]C, H,-1.j,1,. .lC, [-1 JJ.-j.- .]C and D. -- ]£)• 
45 straining the phase response of the equalizer B103 to an integer multiple of n at the preamble frequency (1/4T) implies 
that the imaginary component of its frequency response is zero. 



so 



55 



0 

-1 

0 

1 



^ C • = 0 



If the imaginary component of the frequency response is constrained to zero, as described above, then constraining the 
magnitude of the equalizer to g at the preamble frequency (1/4T) implies that the real component of the frequency 
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response equals g. 



1 

0 
-1 
0 



= C -v^ = g 



Therefore, the equalizer's coefficients must be constrained to satisfy the following two conditions: 

C'^-V^ =Oand C^-Vg = g. 

IS 

The above constraints are achieved by multiplying the computed gradient • e^ by an orthogonal projection operation 
Pv^vg ^ as part of a modified LMS algorithm C1 14. 

To understand the operation of the orthogonal projection operation, consider an equalizer that comprises only two 
coefficients: Cq and as shown in FIG. 10. The phase constraint condition • V ^ =0 implies that the filter coeffi- 

20 cient vector must be orthogonal to V-, . When using an unmodified LMS algoTithm to update the filter coefficients, the 
orthogonal constraint is not always satisfied as shown in FIG. 10. The present invention, however, constrains the fitter 
coefficients to a subspace <C ) which is orthogonal to Vi by multiplying the gradient values X|^ • eK by a projection oper- 
ation Pv^-^. where the null space of the projection operation Pv/ is orthogonal to (C ). The updated coefficients corre- 
spond to a point on the orthogonal subspace <C) closest to the coefficients derived from the unmodified LMS algorithm 

25 as shown in FIG. 10. 

Similar to the phase constraint projection operation Pv-,-*-, a second orthogonal projection operation Pv2^ constrains 
the filter coefficients such that the coefficient vector satisfies the above gain constraint: C ^ • V g = g • The combined 
orthogonal projection operation Pv^vg-^ eliminates two degrees of freedom in an N-dimensional subspace where N is 
the number of filter coefficients (i.e., the orthogonal projection operation PviV2^ has a rank of N-2). 
30 An orthogonal projection operation for and V2 can be computed according to 

Pv/ = l-Pv, = l-V,(V/v,)-\/ (20) 

where Pv ^Vg ^=Pv , ^ * PV2 since Vi is orthogonal to Vg- The orthogonal projection operation PviV2^ computed 
35 using the above equation for an equalizer comprising ten filter coefficients is a matrix 
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4 



50 



The above matrix Pv^V2^ is an orthogonal projection matrix scaled by 5 (multiplied by 5) so that it contains integer 
valued elements which simplifies multiplying by • e^ in the LMS update equation, 

W,,, = W,-M-PviV2^-(X,.eK)- (21) 

The scaling factor is taken into account in the selection of the gain value m Constraining the gain to g and the phase to 
kit at the normalized frequency of 0.5 simplifies implementing the orthogonal projection matrix Pv^ V2-'-: half of the ele- 
ments are zero and the other half are either +1.-1, or +4. Thus, multiplying the projection matrix Pv^ V2-^ by the gradient 
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values X|^ * e^^ requires only shift registers and adders. 

The ACQ/TRK signal shown in FIG. 8B disables adaptation of the FIR filter during acquisition, that is. while acquir- 
ing the acquisition preamble 68 shown in FIG, 2B. Thus, the adaptive equalizer B103 adapts only after acquiring the 
acquisition preamble 68. 

Reduced Cost Orthogonal Constraint Matrix 

Even though the above orthogonal projection matrix Pv^ V2-^ has a simple structure wherein half of the elements are 
zero, it may not be cost effective to directly implement it due to the significant number of shift and accumulate operations 
necessary to compute Pv ^ v g ^ • (X • e . In order to reduce the cost and complexity, an alternative embodiment of 
the present invention decimates the modified LMS adaptation algorithm C114 as illustrated in FIG.s 11 A. 11 B, 110 and 
11D. 

A mathematical basis for the circuit of FIG. 1 1 A will be discussed before describing the details of operation. Refer- 
ring again to the above equation (20), 

Pv/ = I-Pv,=l-V,(v/V,)-^V/ 



by combining the above Vf and V2 vectors into a Nx2 matrix 



20 



25 



V = 



+1 


0 


0 


-1 


-1 


0 


0 


+1 



(those skilled in the art will recognize that shifting the time base provides four alternatives for the V matrix) then the 
30 operation (V^^Vx)"^ of equation (20) reduces to 



35 



1 0 
5 1 
0 5 



■I' 



Thus, equation (20) reduces to 



Pv-^= l-l/5-Vv"^. 



40 



Multiplying both sides of equation (22) by 5 provides 

5-Pv-^=5-l- W"^. 

45 Referring again to equation (21), 

Wk.i = W,-M-PviV2^.(X,-e,) 
setting Xj^-e^ = gj^ reduces equation (21) to 

Defining 8 = v^ • g^ 



(22) 



(23) 



50 



(24) 
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+10-10 
0-10+1 



70 



go 

g3 



5, 



(25) 



IS 



25 



that is. 8o=go-g2+ -- a^K* 61=-g1+g3- 
then computing • 





+1 0 




+5o 




0 -1 




So" 




-5, 


20 


-1 0 








-So 




0 +1 



















and computing 5-g^-w^-g^ provides 



(26) 



30 


+5o 




5go-5o! 




-5i 








-So 




5g2+S^ 


35 


+6, 




5g3-5, 



(27) 



40 



FIG. 1 1 A implements the above equations (25). (26). (27) and ultimately equation (21) in order to update the coef- 
ficients of the equalizer filter CI 00 according to the modified LMS adaptation algorithm of the present invention. To fur- 
ther reduce the implementation cost, the circuit of FIG. 1 1 A decimates the adaptation algorithm by the number of filter 
coefficients (10 in the example shown); that is, the adaptation algorithm operates only on every tenth sample value 
45 rather than on every sample value and updates only one filter coefficient per clock period. This is illustrated by expand- 
ing equation (24) 



50 



55 



w 

^ (lc+ 1)0 




w 




W 




W 




W 

^(k + 3)2 




^(k + 2)2 




w 




w 

^ (k + 3)3 





5g(k + o)o-S(^^o)o 

5g(k + 3)3-S(k.^.3), 



(28) 



where k+i is the clock period and g = X ^ • e j . Thus, at each clock period, the next gradient value g^jjj can be com- 
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puted by multiplying a sample value latched every tenth clock cycle by the next error value ej. The new gradient g^jjj is 
then used to update the corresponding filter coefficient \N^^ of equation (28). 

Referring now to the circuit shown in FIG. 1 1 A, the sample values X^^ are input into the discrete time equalizer filter 
ClOO and the equalized sample values input into the interpolator B122 similar to FIG. 8B. A serial-to-parallel circuit 

5 C140 converts the interpolated sample values into even and odd subsequences Y2k+x ancJ Y(2k.i)+T. where the 
notation 2k indicates that two new interpolated sample values are output from the serial-to-parallel circuit C140 at every 
other sample period. A slicer CI 42 generates the corresponding even and odd estimated subsequences ~Y2k+x ^^d 
-Y(2k-i)+T which are subtracted from the interpolated sample values at respective adders (C144^,C1442) to form an 
even and odd sample error sequences e2k+T anci e(2k-i)+T- An error value interpolation circuit C146, similar to that of 

70 FIG. 8B, generates the even and odd sanple error sequences e2k and e(2k-i) which are synchronized to the A/D sample 
rate. 

As mentioned above, the circuit of FIG. 1 1 A decimates the adaptation algorithm of the present invention by the 
number of coefficients in the equalizer filter C100. For the example shown, the equalizer filter CI 00 has 10 filter coeffi- 
cients; accordingly, a decimate by 10 circuit C148 loads a sample value into a delay register CI 50 every tenth sample 

15 period. Thereafter, the output of delay register CI 50 is represented by X2k-2j-i where; j=1 ->5 incremented by 2 at every 
other sample period. The output of the delay register CI 50 X2k-2j-i 's multiplied by the sample errors e2k and e(2k-i) a* 
respective multipliers (C152-|,C1522) to form the gradient values g2j and g2j-i used in equation (24). 

The gradient values g2j and g2j-i are then shifted into respective shift registers (CI 54^, CI 542). To implement equa- 
tion (27), the gradient values g2j and g2j-i are multiplied by respective alternating ±1 (C156t.C1562) and accumulated 

20 in respective accumulators (CI 58^, CI 583). After accumulating 5 gradient values in each accumulator (C158i»C1582)« 
the outputs of the accumulators (which represent Sq and of equation (25)) are latched by respective decimate by 5 
circuits (C160i,C1602) and the accumulators (CI 58^, CI 582) are reset. The values 60 and 61 are then multiplied by 
respective alternating ±1 (C162^.C1622) to implement equation (26). The gradient values g2j and g2j-i at the outputs of 
the shift register (C154^,C1542) are multiplied by 5 (C164^,C1642) and the values 80 and 6^ are subtracted therefrom 

25 at adders (CI 66^. CI 662) in order to implement equation (27). 

To finish the adaptive update algorithm (i.e., to implement equation (28)). the output of adders (C1 66.| .CI 662) are 
scaled by a gain factor ^ (CI 68^, CI 682) which is reduced by a factor of 5 to account for the scaled up projection oper- 
ator. The output of the gain factor ^ (CI 68i.C1682) is subtracted at adders (C170i .CI 7O2) from the corresponding filter 
coefficient (W2j,i.W2^ selected by a multiplexor (C172i.C1722) from a bank of registers (C174i,C1742). The ADAPT2j 

30 signal selects the 2j*" coefficient from the bank of registers (C174^,C1742) for updating. After subtracting the update 
value, the updated filter coefficient (C176i,C1762) is restored to the bank of registers (C174^.C1742) and used by the 
equalizer filter C100 during the next clock period to equalize the sample values according to its updated spectrum. 

Decimating the update algorithm by the number of filter coefficients, as in FIG. 1 1 A, decreases the implementation 
complexity but at the cost of slowing convergence of the equalizer filter C100 toward an optimum setting. This is 

35 because the decimated update algorithm does not use all of the sample values to compute the gradients gj. FIGs. 1 1 B 
and 110 show an alternative embodiment of the decimated update algorithm of the present invention wherein more of 
the sample values are used to compute the gradients gj. which adds complexity but improves performance because the 
equalizer filter CI 00 converges faster. 

The circuit in FIG. 118 operates similar to the circuit in FIG. 1 1 A described above except for the addition of respec- 

40 tive gradient averaging circuits (CI 76^, CI 762) which compute an averaged gradient g2j-i and g2j over several sannple 
values. 

^ N - 1 

9 2kj ~ N 2 ®2k-n * ^ 2k-2j-n (29) 
"^5 n = 0 



where N is a predetermined number of sample values out of the number of equalizer filter taps. In one embodiment, all 
of the sample values could be used (i.e., a number of sample values equal to the nunnber filter taps), or in an alternative 

50 embodiment, a decimated number of sample values could be used in order to reduce the complexity and cost. 

FIG. 1 10 shows an embodiment of the gradient averaging circuits (CI 76^.01762) wherein the number of sample 
values is decimated by 2; that is, 5 out of the 1 0 sample values corresponding to the 1 0 equalizer filter taps are used to 
connpute the averaged gradients gaj-i and g2j. 

In operation, a decimate by 2 circuit CI 78 of FIG. 1 1 B stores every other sample value in the delay register CI 50. 

55 The output of the delay register Cl 50 is multiplied (0180^-1805) by the error value e2k and delayed versions of the error 
value e2k.n. and the results are accumulated in accumulators (C1 82^-1825). After accumulating five gradients, the con- 
tents of the accumulators (C182i-1825) are transferred to registers (0184^-1845). and the accumulators (C182^-1825) 
are cleared. Then, at every other sample period, the contents of registers (0184^-1845) are shifted from left to right (i.e.. 
01845 = OI844; 01844 = OI843: etc.) and the output of register OI845 is the averaged gradient g2j output by the gra- 
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10 



dient averaging circuit (C176i,C1762). The averaged gradients gsj and ggj.! (C186i,Cl862) are then used to update the 
coefficients of the equalizer filter C100 using circuitry CI 88 in the same manner as described with reference to FIG 
11 A. 

FIG. 1 1 D illustrates yet another embodiment of the present invention which further reduces the cost and complex- 
ity, as conrpared to FIG. 1 1 A. by updating the even and odd filter coefficients sequentially. That is. during the first N peri- 
ods (where N is the number of filter coefficients) the circuit of FIG. 11 D computes the coefficient update values of 
equation (28) for the even filter coefficients (Wq, W2. W4...). and then during tiie next N sample periods it computes the 
update values for the odd filter coefficients (W^. W3, W5,...). A decimation circuit CIdO decimates the error value e^ by 
two. and approximately half the circuitry as that of FIG. 1 1 A is used to compute the update values. The decimate by two 
circuit CI 90 is actually syncopated; that is. it outputs the error values for k=0, k=2. k=4, k=6. k=8. and then outputs the 
error values for k=9. k=i 1 . k=l3. k=15. k=l7 (assuming the adaptive filter CI 00 comprises ten filter taps). 

Mathematically, operation of the update circuit in FIG. 1 1D can be described by, 



75 



20 



25 



30 



w 

^ (k+ 1)0 




w 

^ (k + 0)0 


w 




^(k + 2)2 


W 




W 

^ (k + 4)4 


w 




w 

^ (k + 6)6 


^(k + 9)8 




^(k + 8)8 


w 

(k+ 10) I 




w 

^ (k+9)l 


W 

^ (k+ 12)3 




W 

^ (k-K 11)3 


W 

^ (k+ 14)5 




W 

^ (k+ 13)5 


W 

^ (k+ 16)7 




W 


W 

(k+ 18)9 




W 

(k + 17)9 



^S(k + 0)O"S(k + 0)O 



5g 



(k + 2)2 • "(^ + 2)0 
5g(V4,4)4-S^j^^^^^ 

5g(k + 6)6 + S(k + 6)0 
5g(k + 8)8-5(t + 8)0 
5g(k*9)l + S' 



()c + 9) I 



5g(k^.Il)3-S(^^,jjj 



5g 



+ 5. 



(k+ 13)5 ■ ^(k + 13) 1 
5(k*l5)7-5^j^^ 



5g 



_-'&(k+ 17)9 ^^()c+ 17)1^ 



35 assuming xt>e adaptive filter ClOO comprises 10 filter taps. 

Obviously, the embodiment of FIG. 11D will decrease the performance of the adaptive filter ClOO due to the 
decrease in coovef gence speed. However, the gradient averaging circuit of FIG. 1 1 G can be used to improve the per- 
formance 0* the ofcuit in FIG. 1 1D. similar to the embodiment of FIG. 1 1B. Thus in the embodiments of FIG. 1 1 A-11D. 
the periofmarce vefsus cost and complexity varies-the preferred configuration is selected according to the require^ 

40 ments of tt>e usei 

Atthough ihe mtei poiated timing recovery and adaptive equalizer of the present invention have been disclosed in 
relation to a d.O PR4 read channel, the principles disclosed herein are equally applicable to other types of sampled 
amplitude iei*d chdfinefe including d=1 EPR4 or EEPR4 read channels. In a d=1 read channel, for example, ttie slicer 
B141 of rK3 B4A IS leplaced by a pulse detector as described in the above reference U.S. Patent No. 5,359,631. 
45 Furthefmoie the particular constraint frequency of 1/4Ts used in the disclosed embodiment is not intended to be 

limiting Othei constraint frequencies could be used without departing from the scope of the present invention. For 
example, a 3T pfcambio could be used in which case the constraint frequency would be 1/6Ts (if constraining to the 
preamble ficquency) 

Still further, those skilled in the art will appreciate the many obvious modifications that are possible to tiie adaptive 
50 equations disclosed herein. For example, a shift in tiie time base would change the FIR filter's magnitude and phase 
response at the constraint frequency which would result in different V, and mati-ices. Also, the modified U^S algo- 
rithm, 

Wk.i = WK-M-PviV2^.(X,.e,) 

55 

could be implemented after rearranging terms, 

Wk.i = Wk-M-e^-(Pv,v2-'-)<^) 
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or 

WK.i=W,-n-X,-(Pv,v/.e,). 

5 These, and like rrxxJifcations, are within the scope of the present invention. 

The objects of the invention have been fully realized through the embodiments disclosed herein. Those skilled in 
the art will appreciate that the various aspects of the inversion can be achieved through different embodiments without 
departing from the essential function. The particular embodiments disclosed are illustrative arxl not meant to limit the 
scope of the invention as appropriately construed by the following daims. 



Table 1 



Channel 


Transfer Function 


Dipulse Response 


PR4 


{1-D)(1+D) 


0,1.0, -1.0. 0. 0, ... 


EPR4 


(1-D)(1+D)2 


0. 1, 1. -1, -1. 0. 0. ... 


EEPR4 


(1-D)(1+D)3 


0.1,2. 0, -2.-1.0. ... 
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Table B2 



5 



15 



25 



35 



40 



SO 



x-32/Ts 


c(-2) 


c(-2) 


C(0) 


c(1) 


c(2) 


C(3) 


0 


0.0000 


-0.0000 


1-0000 


0.0000 


-0.0000 


0-0000 


1 


0.0090 


O.0231 


0.9965 


0.0337 


-0.0120 


0.0068 


2 


0.0176 


-0.0445 


0.9901 


0.0690 


-0.0241 


0.0135 


3 


0.0258 


-0.0641 


0.9808 


0.1058 


-0.0364 


0.0202 


4 


0.0335 


-0.0819 


0.9686 


0.1438 


-0.0487 


0.0268 


5 


0.0407 


-0.0979 


0.9536 


0.1829 


-0-0608 


0-0331 


6 


0.0473 


-O.1120 


0.9359 


0,2230 


-0.0728 


0.0393 


7 


0-0533 


-0.1243 


0.9155 


0.2638 


-0.0844 


0.0451 


8 


0.0587 


-0.1348 


0.8926 


0.3052 


-0.0957 


0.0506 


9 


0.0634 


-0.1434 


0.8674 


0.3471 


-0.1063 


0.0556 


10 


0.0674 


-O.1503 


0.8398 


0.3891 


-0.1164 


0.0603 


11 


0.0707 


-0.1555 


0.8101 


0.4311 


-0.1257 


0.0644 


12 


0.0732 


-0.1589 


0.7784 


0-4730 


-0.1341 


0.0680 


13 


0.0751 


-0.1608 


0.7448 


0-5145 


-0.1415 


0.0710 


14 


0.0761 


-0.1611 


0.7096 


0.5554 


-0.1480 


0.0734 


15 


0.0765 


-0.1598 


0-6728 


0.5956 


-0.1532 


0.0751 


16 


0.0761 


-0.1572 


0.6348 


0.6348 


-0.1572 


0.0761 


17 


0.0751 


-0.1532 


0.5956 


0.6728 


-0.1598 


0.0765 


18 


0.0734 


-0.1480 


0.5554 


0.7096 


-0.1611 


0.0761 


19 


0:0710 


-0.1415 


0.5145 


0.7448 


-0.1608 


0.0751 


20 


0.0680 


-0.1341 


0.4730 


0.7784 


-0.1589 


0.0732 


21 


0.0644 


-0.1257 


0.4311 


0.8101 


-0.1555 


0.0707 


22 


0.0603 


-0.1164 


0.3891 


0.8398 


-0.1503 


0.0674 


23 


0.0556 


-0.1063 


0.3471 


0.8674 


-0.1434 


0.0634 


24 


0.0506 


■0.0957 


0.3052 


0.8926 


-0.1348 


0.0587 


25 


0.0451 


-0.0844 


0.2638 


0.9155 


-0.1243 


0.0533 


26 


0.0393 


-0.0728 


0.2230 


0-9359 


-0.1120 


0.0473 


27 


0.0331 


-0.0608 


0.1829 


0.9536 


-0.0979 


0.0407 


28 


0.0268 


•0.0487 


0-1438 


0.9686 


-0.0819 


0.0335 


29 


0.0202 


-0.0364 


0.1058 


0-9808 


-0.0641 


0.0258 


30 


0.0135 


-0.0241 


0.0690 


0.9901 


-0.0445 


0.0176 


31 


0.0068 


-0.0120 


0.0337 


0.9965 


-0.0231 


0.0090 
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A sampled amplitude read channel for reading digital data from a sequence of discrete time sample values gener- 
ated by sampling an analog read signal from a read head positioned over a magnetic medium, comprising: 
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(a) a sampling device for sampling the analog read signal to generate the discrete time sample values; 

(b) an adaptive equalizer, responsive to the discrete time sanple values, for generating equalized sannple val- 
ues according to a target response: 

(c) interpolated timing recovery for generating interpolated sannple values; and 

5 (d) a discrete time sequence detector for detecting the digital data from the interpolated sample values. 

2. The sampled amplitude read channel as recited in claim 1 . wherein: 

(a) the adaptive equalizer comprises an FIR filter comprising a plurality of filter coefficients; 

10 (b) the adaptive equalizer operates according to a least mean square algorithm responsive to an eirror value ek 

computed as a function of an output of the FIR filter and an estimated ideal value; 

(c) the adaptive equalizer comprises an interpolation circuit, responsive to the interpolated sample values, for 
generating the error value e|^ synchronous with the discrete time sample values; and 

(d) the least mean square algorithm uses the error value to ac^ust the filter coefficients of the FIR filter. 

15 

3. The sampled anplitude read channel as recited in claim 2. wherein the error value e^^ Is computed relative to a dif- 
ference t>etween an output of the interpolated timing recovery and a sample value estimator responsive to the inter- 
polated sample values. 

so 4. The sampled arrplrtude read channel as recited in claim 2, wherein the error value e^ is computed relative to a dif- 
ference between an output o1 the FIR filter and an output of the interpolation circuit. 

5. The sampled annplitude read channel as recited in claim 2. wherein: 
25 (a) the least mean square algorithm is: 

WK.i=WK-M-X,-e,: 

(b) Wi. are the fitter coefficients of the FIR fflter; 
30 (c) ft IS a predetermined gain; and 

(d) IS a discrete time sample value. 

6. The sampled amplitude read channel as recited in claim 2. wherein the adaptive equalizer further comprises a con- 
straint circuit lof constraining a frequency response of the adaptive equalizer at a predetermined constraint fre- 

35 quency 

7. The sarrpied arrplrtude read channel as recited in claim 6, wherein the constraint circuit comprises an orthogonal 
projection opef ation for constraining the filter coefficients relative to the predetermined constraint frequency. 

40 8. The sampled amplitude read channel as recited in claim 7. wherein: 

(a) the constraint orcuit operates according to a modified least mean square algorithm, 

Wk.i =W,.M.Pv,V2^-(X,.eJ 



45 
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(b) W*. are the l»ttef coefficients; 

(c) M 10 a prcdctcf mined gain; 

(d) Pv , v^. ic the orthogonal projection matrix; and 

(e) IS a di&crete time sample value. _ 

9. The sampled amplrtude read channel as recited in claim 1 , wherein the interpolated timing recovery further gener- 
ates a data ctoc*< tor clocking the discrete time detector 

10. The sampled amplitude read channel as recited in claim 1 , wherein: 

(a) the analog read signal conrprises pulses modulated by the digital data at a predetermined baud rate; arKi 

(b) the sampling device samples the analog read signal at a rate above the baud rate. 

11. The sampled annplitude read channel as recited in claim 1 , wherein the interpolated timing recovery comprises an 
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interpolator responsive to an interpolation interval x proportional to a time difference between a channel sample 
value and a desired interpolated sample value. 

12. The sampled amplitude read channel as recited in claim 1 , wherein the analog read signal comprises pulses mod- 
ulated by the digital data at a predetermined baud rate, the interpolated timing recovery further conprising; 

(a) a frequency offset generator for generating a frequency offset proportional to a frequency difference 
between a sampling clock and the baud rate; and 

(b) a mod-Ts accumulator for accumulating, modulo-Ts, the frequency offset A/ to generate an interpolation 
interval x where Ts is a predetermined sample period of the sanrpiing clock. 

13. The sampled amplitude read channel as recited in claim 12. wherein the frequency offset generator comprises: 

(a) a phase error detector for detecting a phase error A0 between an interpolated sample value and an esti- 
mated sample value; and 

(b) a loop filter for filtering the phase error AO to generate the frequency offset Af. 

14. The sampled amplitude read channel as recited in claim 11. wherein the interpolator comprises an actual interpo- 
lation filter having an actual impulse response h{k) which approximates an ideal impulse response of an ideal Inter- 
polation filter: 

sinc(7t • (k-x/T s)); where: 

k is a time index; and 

Ts is a sample period of the sampling clock. 

5. The sampled amplitude read channel as recited in claim 11. wherein the interpolator comprises an actual interpo- 
lation filter having an actual impulse response h(k) generated by minimizing a mean squared error between a fre- 
quency response of the actual interpolation filter and an ideal frequency response corresponding to an ideal 
interpolation filter. 

6. The sampled amplitude read channel as recited in claim 1 1 , wherein: 

(a) the interpolator comprises an actual interpolation filter; and 

(b) a plurality of coefficients k oi the actual interpolation filter are computed in real time as a function of x. 

7. A method for reading digital data from a sequence of discrete time sample values generated by sampling and ana- 
log read signal from a read head positioned over a magnetic medium, comprising: 

(a) sampling the analog read signal to generate the discrete time sample values; 

(b) adaptively equalizing the discrete time sample values to generate equalized sample values according to a 
target response; 

(c) interpolating the equalized sample values to generate interpolated sample values; and 

(d) detecting the digital data from the interpolated sample values. 

8. The method for reading digital data as recited in claim 17, wherein the step of interpolating sample values com- 
prises the step of generating a data clock for clocking the discrete time sequence detector. 
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